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ABSTRACT.  This  report  gives  the  derivation  of  equations  for  the  ray  trace 
analysis  of  an  mtenna  consisting  of  a  parabolic  dish  with  a  phased  array  in 
front  of  it.  The  deri  /alien  gives  the  equations  for  ray  tracing  in  a  form  stri.- 
/  Die  for  digital  computer  use.  A  number  of  rays  from  a  plane  wavefront  can 
be  traced  through  the  antenna  rystein  to  a  feed  point  on  the  antenna  axis.  The 
rays  are  straight  lines,  aid  specular  reflection  is  assumed  at  the  dish.  The 
phased  array  returns  the  reflected  rays  to  the  feed  point  (rvhicb  is  assented  to 
lie  between  the  dish  and  phased  array). 

The  major  features  provided  by  the  analysis  and  the  computer  program  are  the 
relative  phases  and  path  lengths  of  the  rays  from  the  wavefront  to  the  feed 
point  (tabulation  and  computer  plots),  the  total  number  of  rays,  the  number 
obscured  by  the  phased  array,  the  number  that  arrive  at  the  phased -may  plane 
and  lie  within  the  phased  array  itself,  a  computer -gene  rated  plot  of  thest-rtmg 
points  cm  the  wavefront  of  the  rays  used  in  tracing,  a  plot  of  the  locations  of 
those  rays  that  arrive  at  the  phased -array  plane,  and  an  outline  of  the  phased 
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FOREWORD 

Existing  C-band  instrumenution  radar  sets  acquire  and  track  targets  by  m-tchanicaily 
positioning  the  antenna  bat,  because  of  the  large  mechanical  mass  that  must  be  moved, 
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without  having  to  move  the  luge  mechanical  mass  of  the  reflector  and  feed  assembly. 
Phased -array  antennas,  which  provide  electronic  beam  steering,  have  been  demonstrated 
to  be  practical  in  other  applications  and  coaid  be  used  to  advantage  in  instrumentation 
radar  sets.  The  Naval  Weapons  Center  is  investigating  a  system  consisting  of  a  small 
phased  array  placed  in  front  of  a  parabolic  reflector.  This  system  is  expected  to  pro¬ 
vide  an  electronically  steered  beam  at  relatively  low  cost  and  -ith  little  modification 
to  existing  sets.  Before  radar  modifications  are  undertaken,  however,  NWC  is  investi¬ 
gating  the  advantages  and  possible  difficulties  inherent  in  the  proposed  system. 

This  report  derives  the  equations  for  the  ray-trace  analysis  of  such  an  antenna  with  a 
phased  array.  An  oatiine  of  *  digiicl  computer  program  that  performs  the  ray-trace 
analysis  is  included. 
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INTRODUCTION 


The  analysis  given  here  provides  the  equations  for  tracing  rays  from  a 
plane  wavefront  through  an  antenna  system  to  a  specified  feed  point.  The  an¬ 
tenna  system  consists  of  a  parabolic  dish  (paraboloid  of  revolution  about  the 
axis  of  the  parabola),  a  circular  phased  array  perpendicular  to  the  paraboloid 
axis  and  centered  on  thi«  axis,  and  a  feed  point  on  the  axis.  The  rim  of  the 
dish  lies  in  a  plane  perpendicular  to  the  axis  and  is,  therefore,  a  circle  cen¬ 
tered  on  the  axis.  The  feed  point  lies  between  the  dish  and  the  phased  array. 
The  phased  array  is  assumed  to  lie  outside  the  dish. 

A  ray  normal  to  the  wavefront  is  traced  from  the  wavefront  to  the  dish. 

It  is  specularly  reflected  from  the  dish  to  the  phased  array,  which  directs  it 
to  the  feed  point.  Those  rays  which  are  not  reflected  by  the  dish  are  omitted 
and  those  which  are  intercepted  by  the  back  of  the  phased  array  can  be  omitted. 
Further,  those  rays  which  are  not  reflected  to  the  face  of  the  phased  array  can 
be  rejected  or  noted.  Conditions  for  multiple  reflection  are  given  so  that  such 
rays  can  be  rejected  also.  The  total  geometric  path  length  from  the  wavefront 
to  the  feed  point  is  found  and  can  be  related  to  the  relative  phase. 

In  the  process  of  ray  tracing  the  projections  of  the  dish  rim  and  the  phased 
array  rim  on  the  wavefront  plane  and  the  image  of  the  dish  rim  on  the  plane  of 
the  phased  array  are  found.  Plots  of  these  curves  and  the  intersections  of  the 
rays  on  the  phased  array  plane  are  useful  in  the  evaluation  of  the  ray  tracing. 

Since  the  ray  tracing  equations  are  complicated  functions  of  the  constants 
defining  the  antenna  geometry,  it  is  a  practical  necessity  to  use  a  large-scale 
digital  computer  for  this  ray  tracing.  The  development  of  the  equations  in  the 
body  of  the  report  is,  therefore,  oriented  to  numerical  analysis.  A  summary 
of  the  equations  suitable  for  the  digital  computer  is  given  in  Appendix  A,  and 
examples  of  computer  output  are  given  in  Appendix  B. 
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A  computer  program  for  ray  tracing  has  been  written  in  FORTRAN  V  for 
the  UNTVAC  1108  with  a  Stromberg  DatagraphiX,  Inc. ,  4060  plotter.  This 
program  (WAVE)  is  available  from  Commander,  Naval  Weapons  Center  (Code 
3037),  China  Lake.  California  93555. 

Thanks  are  due  John  Anderson  and  Lorene  Paschal,  Code  3037,  and  Edlin 
Patterson,  Code  3045,  for  useful  discussions  during  the  course  of  this  work. 
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DERIVATION  OF  THE  RAY  TRACING  EQUATIONS 


The  geometry  of  the  antenna  system  is  shown  in  Fig.  1  and  2.  The  origin  of  the 
xyz  coordinate  system  is  at  the  vertex  of  the  parabolic  dish  and  the  x  axis  is  along 
the  dish  axis.  The  following  notation  is  used: 

2,2., k  unit  vectors  along  x,  y,  and  z,  respectively 

Q  position  vector  of  a  point  on  the  dijh 

n  unit  (inward)  normal  to  the  dish  at  point 

N  unit  normal  to  wavefront  plane  (from  vertex  to  plane) 

N  distance  from  origin  to  xyz  to  wavefront  plane  along  N 

P  position  vector  of  point  on  wavefront 

R  unit  vector  along  reflected  ray 

J,K  unit  vectors  along  Y  and  Z  coordinate  axes  lying  in  the  wavefront 
plane 

NN  the  position  vector  of  the  origin  of  the  YZ  coordinate  system 
Y,  Z  coordinates  of  point  on  wavefront  in  YZ  coordinate  system 
A  radius  of  rim  of  dish 

F  focal  distance  of  dish 

D  distance  (along  x)  from  origin  to  phased  array 

R0  radius  of  phased  array 

f  distance  (along  x)  from  origin  to  feed  point 

yp,z  position  coordinates  of  intersection  of  reflected  ray  with 
^  phased  array 

Lt  path  length  from  wavefront  to  dish 

Lj  path  length  of  reflected  ray  from  dish  to  phased  array 

Lg  path  length  from  phased  array  to  feed  point 

Other  notation  is  defined  as  it  is  used. 
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FIG.  1.  Antenna  System  Showing  the  Three  Segments,  L^,  L^,  L^,  of  the  Path  From  the  Wavefront 
to  the  Feedpoint. 
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FIG.  2.  Projection  of  the  Incident  Ray  Onto  the  Dish  and  Specular  Reflection. 
Vectors  N,  P-NN,  2j  P  lie  in  a  common  plane. 
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The  equation  of  the  parabolic  dish  is 


where  Qy. 
position  vector 


Q2  +  Q2  =  4FQ 
y  z  x 

and  Q„  are  the  coordinates  of  a  point 

j Q 

9.  =  Q 

\q 


=  i Q  + iQ 

-  x  jL  y 


on  the  dish. 


+  kQ 
—  z 


Therefore,  the 


(1) 


(2) 


The  unit  vector  N,  which  defines  the  orientation  of  the  wavefront  plane  has  com¬ 
ponents  that  are  just  the  direction  cosines  of  the  normal  to  the  plane.  Since  the 
antenna  system  has  rotational  symmetry  about  the  x  axis,  it  is  sufficient  to  keep  N 
in  the  xy  plane  and  define  its  direction  by  a  single  angle  a.  Then 


N 


(3) 


where  a  is  the  angle  from  the  x  axis  to  N.  For  the  present  purpose  0  <  a  <  90°  is 
the  only  range  of  interest.  Reflection  from  the  convex  side  of  the  dish  is  not  required. 
Negative  a  is  not  required  since  the  system  is  mirror  symmetric  about  the  xz  plane. 


Projection  of  Dish  and  Phased  Array  Rims  on  Wavefront 

With  a  plane  wavefront,  the  rays  from  it  to  the  dish  are  all  parallel  to  N.  The 
ray  tracing  from  the  wavefront  to  the  dish  consists  of  finding  the  position  of  the  ray 
which  starts  with  coordinates  Y  and  Z  on  the  wavefront.  Only  the  region  on  the  wave- 
front  within  the  projection  of  the  rim  of  the  dish  need  be  considered.  The  projection 
of  the  rim  is  made  parallel  to  N.  The  projection  of  the  rim  of  the  phased  array  can 
obscure  part  of  this  area;  rays  within  the  intersection  of  the  two  projections  can  either 
be  omitted  or  noted.  It  is  of  interest  to  trace  such  rays  to  show  the  extent  of  the 
"shadow". 


The  computational  procedure  is  to  trace  rays  from  a  rectangular  grid  of  points 
on  the  wavefront.  Only  the  points  which  lie  in  or  on  the  projection  of  the  dish  rim 
are  used.  It  is  now  useful  to  find  the  bounds  of  this  region  on  the  wavefront. 
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(8) 


Z  =  a  -  K  =  Qz 


anu  the  equation  of  the  rim  of  the  dish  is 


Q*y  +  Q*z  =  4FQx  =  A2  (2) 

Combining  Eq.  (7),  (8),  and  (9)  gives 

Z2  +  +  /cos2  a  =  A2  (10) 

This  equation  is  the  projection  of  the  rim  on  the  wavefront  plane  and  is  an  ellipse  with 
center  at  Z  =  0,  Y=  -  A2  sina/4F.  The  semimajor  axis  is  A,  and  the  semiminor  axis 
is  Acosa.  This  is  illustrated  in  Fig.  3. 
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For  computer  plots,  it  is  convenient  to  put  ^q.  (10)  in  parametric  form. 

At  the  rim. 


Q  =  A  sin  0 

'j 

Q  =  A  cos  0 
z 

Then  Eq.  (7)  and  (8)  become 


Y 


A  sin  6  coso 


A* 

4F 


sin  a 


Z  -  a  cos  a 


(10a) 


\ 

i 


i 

I 


for  the  projection  of  the  dish  rim  an  the  wavefront  plane.  It  is  sufficient  to  use 
-a/2  <fl  <  a/2  and  use  symmetry  about  the  Y  axis  to  plot  for  Z  <  0. 

Next,  the  projection  of  the  rim  of  the  phased  array  is  found.  Since  the  phased 
array  is  circular,  the  equations  of  the  rim  are 


x  =  D  and  y2  +  z2  =  r2 

p  p  P 


fll) 


and  the  projection  of  the  rim  position  vector  Dj_  +  yj_+zk=r  onto  the  wave- 
front  plane  gives  ^  ** 


Y  =  r  ■  J  =  -  D  sin  a  +  y  cos  a 
-  P 


Z  =  r  K  = 


z 

P 


(12) 


Combining  Eq.  (11)  and  (12)  gives  the  equation  of  the  projection  of  the  phased  array 
rim, 


Z2  +  (Y  +  Dsinaf/cos2a  =  R§ 


(13) 


This  is  an  ellipse  with  the  center  at  Z  =  0,  Y  =  -  D  sina,  and  the  semimajor  and 
semiminor  axes  are  R3  and  Rg  cos  a,  respectively.  See  Fig.  3. 
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Again,  for  iMhn  computer  plots  it  is  convenient  to  have  Eq.  (13)  in  parametric 
form.  Following  the  scheme  to  get  Eq.  (10a) 

\  =  R|  rtiil  0  COS U  -  u  SSalCt 

Z  =  R*  cos  8  (13a) 

-  v  <  S  <  f 

for  the  coordinates  of  the  projection  of  die  phased  array  rim. 


Points  on  Wavefront  to 


For  computational  purposes,  a  grid  of  points  on  the  YZ  plane  is  used  to  start  the 
ray  tracing.  Only  those  points  on  or  within  the  ellipse  of  Eq.  (10)  are  used. 

If  rays  that  start  on  points  along  a  line  parallel  to  the  T  axis  are  required,  the 
limits  cm  Y  imposed  by  Eq.  (10)  are 


(A2  -Z2)^  cos  a  -  A^pl~  ^  Y  £  (A2 -Z2)1''2  cos  a  - 


A*  sina 
4F 


where 


Z2  ^  A1 


Considering  Eq.  (13),  if 


R?  :£  Z2  <  A1 


Then  the  full  range  of  Y  given  in  (14)  is  used.  On  the  other  hand,  for 


Z2  <  R? 


all  Y  in  the  range 
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<R2-zV^  cosa  -  D  sin  a  <  Y  <  (Rl-Z)1^  cos  a  -  D  sin  a 


(18) 


are  within  the  projection  of  too  phased  array. 

Since  the  antenna  system  is  mirror  symmetric  about  the  xy  plane  and  N  is  in 
the  xy  plane,  it  is  sufficient  to  use  Z  £  0  only. 

For  programming  the  computer  it  is  convenient  to  have  the  coordinates  of  points 
on  the  wavefront  plane  indexed  by  integers  anH  to  have  the  limits  in  terms  of  these 
indices.  For  the  Z  coordinate  in  the  range 


0  <  Z  <  A 


Z(J)  =  A  Z  *  (J-l),  J  =  1,2,. ...NZ 


(19) 


(20) 


where  A  Z  is  the  increment  an  Z,  and 


NZ  =  INT 


(A) 


+  1 


(21) 


is  the  number  of  points  on  Z  on  or  within  the  limits  (18)  and  INT  means  "the  integer 
part  of'. 

The  points  on  the  Y  axis  have  the  limits  (14)  (see  Fig.  3).  For  a  given  Z(.T),  the 
number  of  Y  points  from  the  major  axis  to  the  ellipse  with  the  increment  AY  is 


NY(J)  =  INT 


( 


[A2-Z2  (J)]1/2 


AY 


cos 


“) 


(22) 


This  equation  does  not  include  the  point  on  the  major  axis.  The  number  of  points 
along  the  Y  axis  is 


NY(1)  =  INT 


( 


A  cosa  \ 
AY  ) 


(23) 


For  Z  =  0  the  Y  (I)  are 


•f 


& 

! 

i, 
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YO)  =  AY  -  (I-NY(l)-l) 


A2  sing 
4F 


(24) 


where 


I  =  1,2, .. . ,  2NY(1)  +  1 


(25) 


For  other  Z(J),  I  has  the  range 

1  =  1+  NY(1)  -  NY(J), . . .  ,1  +  NY(1)  +  NY(J)  (26) 

Equations  (20)  and  (21)  provide  the  Z  coordinates,  and  Eq.  (22)  through  (26)  provide 
the  Y  coordinates  and  range  of  the  index  I  for  the  Z  that  cover  the  right  half  of  the 
projection  of  the  dish  rim. 

Those  points  Z(J),  Y(I)  lying  within  the  projection  of  the  phased  array  are  now 
readily  found  from  conditions  (16)  through  (18).  Although  the  rays  which  fall  within 
the  phased  array  projection  cannot  reach  the  dish,  it  is  still  useful  to  trace  them  to 
find  the  number  in  this  "shadow  area"  and  to  find  the  image  of  this  shadow  on  the 
phased  array  plane. 


Ray  From  Wavefront  to  Dish 

Next,  the  p>oint  <£  on  the  dish  is  found  given  the  coordinates  X  and  Y  of  a  point 

on  the  wavefront.  Solving  Eq.  (1),  (7),  and  (8)  for  Q  ,  Q  ,  and  Q  gives 

x  y  z 


Q 


X 


Qz+Z2 

y _ 

4F 


2F 


cos  a 
sina 


Z2  sina  +  4FY 
4F2cos2a 


Q 


z 


=  Z 


(27) 
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£ 

p 


' 

i 

jjj 


for  a^O.  The  choice  of  sign  in  the  Qy  equation  is  made  as  follows.  For  Y  =  Z  =  0, 

Q  =  0,  since  the  origin  of  the  YZ  coordinate  system  is  taken  as  the  position  vector 
NN  starting  at  the  origin  of  xyz.  At  the  origin  Qx=Qy=Qz  =  0.  The  negative  sign  in 
th3  Qy  equation  is  therefore  required.  Since  the  antenna  system  has  mirror  sym¬ 
metry  about  the  xz  plane,  only  positive  a  need  be  considered.  This  can  be  seen  from 
thfe  Qy  equation.  If  both  X  and  Y  change  sign,  Qy  changes  sign.  and  Qz  remain  un¬ 
changed.  These  changes  are  just  the  expected  mirror  symmetry. 

For  a=0  Eq.  (7)  becomes  Y=Qy,  then  the  components  of  are 


Z2+Y2 

4F 


Q 


z 


Z 


{28) 


For  8 mall  but  nonzero  a  equation  (27)  are  valid,  but  the  Qy  equations  may  lead  to 
computational  difficulties.  This  equation  can  be  expanded  in  power  series  in  sin  a, 
and  to  second  order 


Q  =  Y/cosa  +  — -  [Z2  +  (Y/cosoqf  ]  Isina  +  ~  Y  2 — sin2a  >  (29) 

y  4Fcosa  '  ' 1  I  2Fcosza  ) 

Using  8  figure  computations  and  typical  numbers  for  Y,  Z,  and  F,  the  Q  equation  of 
(27)  gives  about  5-figure  accuracy  for  a=0.  01  radian.  For  angles  smaller  than 
0.01  radian  Eq.  (29)  will,  in  general,  give  equal  or  better  accuracy. 

Considering  Fig.  2,  and  recalling  Eq.  (3),  the  path  length  from  the  wavefront  to 
point  on  the  dish  is 


Lj  =  N  ~  N  =  N  -  Q  cos  a  -  Q  sin  a 

x  y 


(30) 


Later,  these  path  lengths  will  be  converted  into  relative  phases,  and  N  =  0  can  be 
used,  since  path  difference  is  all  that  is  necessary.  Therefore, 
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L,  =  -  (Q^cosa+Q  sin  a) 


(31) 


In  this  case  the  sign  of  Lj  must  be  retained. 


Reflection  From  Dish  to  Phased  Array 

The  next  path  required  is  that  of  the  reflected  ray  from  the  dish  to  the  phased 
array.  For  this  purpose,  the  unit  normal  n  to  the  dish  surface  is  required.  The 
gradient  of  a  surface  V =V(x,  y,  z)  is  normal  to  that  surface  so,  in  the  case  of  the 
parabolic  dish  V =Q2  +Q2  -4FQx  and  the  gradient 


VV  =  -  4Fi_  +  2Qyl  +  2Qzk 

The  two  unit  normals  are  then 

_  _ VV _ 

"  ±[(4F)2  +  (2Qy)2  +  (2Qz)2]1/Z 

For  the  normal  on  the  concave  side  of  the  dish  the  negative  sign  on  the  square  root  is 
taken.  In  terms  of  components 


The  conditions  for  specular  reflection  are  that  the  angle  of  incidence  and  reflec¬ 
tion  are  equal  and  that  the  incident  and  reflected  rays  and  the  normal  to  the  reflecting 
surface  element  are  coplanar.  In  vector  form  (Fig.  2)  these  conditions  are 


N  •  n  =  n  •  R 


Nx  n  =  n  x  R 


(33) 
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where  R  is  a  unit  vector  in  the  direction  of  the  reflected  r?  y.  Recall  that  all  the  rays 
projected  onto  the  dish  are  parallel  to  N.  To  solve  Eq.  (25)  for  R,  take  the  dot  product 
equation  and  the  y  and  z  component  equations  of  the  cross  product.  In  matrix  form, 
these  are 


1-  n 


fcosa-  n  +  sma-  n 
x  y 


-  cos  a  •  n 


(34) 


cosa  *  n  -  sma  •  n 
y  xy 


the  determinate  of  the  left  square  matrix  is  n^  therefore,  Eq.  (34)  are  not  singular 
unless  nx  =  0.  This  certainly  will  not  occur  for  finite  sized  dishes,  Eq.  (32).  The 
solutions  of  (34)  are 


ft2c?  -1)  cosa  +  2n  n  sin  a' 
'  x  '  x  y 


R  = 


(2n  -i)  sina  +  2n  n  cosa 
y  x  y 


2n  (n  cosa  +  n  sina) 
z  x  y 


(35) 


It  can  be  verified  that  R  is,  indeed,  the  solution  to  the  four  equations  (33)  and  that  it 
is  a  unit  vector. 


Under  some  conditions,  the  reflected  ray  will  not  leave  the  dish  without  a 
second  reflection.  A  necessary,  but  not  sufficient,  condition  is  that 


R  >  0 
x 


(36) 


for  the  ray  to  be  reflected  orly  once.  Reflected  rays  not  meeting  the  condition  (36) 
can  be  rejected  at  this  point.  In  ue  following,  Rx  >0  is  assumed. 

The  equation  of  the  line  with  a  direction  specified  by  R  and  passing  through  the 
point  Q  on  the  dish  is 


rs 
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the  equation  of  the  plane  of  the  rim  of  the  dish  is 


A2  =  4Fx 


(38) 


The  coordinates  of  the  intersection  of  the  line  (37)  and  the  plane  (38)  are,  therefore. 


= 


R 

S-  /Ai 


d  R 


(A2/4F  -  Qx)  +  Q 


R 


(39i 


z.  = 


d  R 


(A2/4F  -Qx)  +  Qz 


The  condition  that  the  reflected  ray  "escape"  from  the  dish  without  a  second  reflec¬ 
tion  is 


y2.  +  z2  <  A2 

d  d 


(40) 


Rays  not  meeting  condition  (40)  are  rejected.  The  conditions,  (36)  and  (40),  are  lor 
the  phased  array  outside  of  the  dish. 


The  equation  of  the  plane  of  the  phased  array  is 


x  =  D 


and  the  intersection  of  the  phased-array  plane  and  the  reflected  ray,  Eq.  (37)  has  co¬ 
ordinates 


R 


y  =  ^  (d-qj  +  Q, 


p  R  '  'x' 

x 


(41) 


z  -  —  (D-Q  )  +  Q 

p  R  x'  z 

x 
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The  condition  for  this  point  to  lie  ou  or  within  the  rim  of  the  phased  array  is 


y2  +  z2  <  R^ 
P  P 


(42) 


Rays  not  meeting  condition  (42)  are  outside  the  phased  array  circle,  but  they  will 
intersect  the  plane  of  the  phased  array. 


Image  of  Dish  Rim  on  Plane  of  Phased  Array 

The  image  of  the  rim  of  the  dish  on  the  phased  array  plane  can  be  useful  in 
evaluating  the  ray  tracing.  This  image  is  the  projection  of  the  rim  along  the  rays 
reflected  from  the  rim.  The  parametric  form  of  the  image  curve  is  found  as 
follows:  At  the  rim. 


Q 

j 

Q. 


Q„ 


A2/4F 

y  . 
rim 

z  . 

rim 


A  sin  9 
A  cos  9 


(43) 


Where  the  parameter  9  has  the  range  -  ir<  9  <  w  to  cover  the  entire  rim.  Substitution 
of  Eq.  (43)  into  Eq.  (32)  and  (35)  gives 


(2^.  -4F2-A2)  sina -4Fy  .  cosa 

rim  rim 


(4F2-A2)cosa  -  4Fy  .  sina 
'  '  Jrim 


2(y  .  sina  -  2Fcosa)  z  . 
rim  _ '  rim 


(4F2  -A2)  cosa  -  4Fy^m  sina 


(44) 


Using  Eq.  (41)  the  coordinates  of  an  image  point  on  the  phased  array  plane  are 
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=  (d-^)  +  y 

*p  \R  I  .  \  4F/  Jrim 

r  \  x/ rim  ' 

/M  |D  ai\  + 

Zp  \Kx)rim  ^  4F'  ‘ri“ 

Rx>«  is  the  necessary  and  sufficient  condition  for  a  ray  reflected  at  the  rim  to  avoid 
multiple  reflections,  and,  if  the  phased  array  is  outside  the  dish,  Rx  >0  insures  the 
ray  intersects  the  plane  of  the  phased  array.  The  denominator  in  Eq.  (44)  comes  from 
Rx  and  must  be  positive  for  Rx  >0.  Thus, 


4F2  -  A2 

y  <  - 

rim  4Ftana 


(46) 


for  0  <  a  <  90°.  For  the  entire  rim  to  have  an  image  in  a  finite  region  of  the  plane 
the  greatest  y  .  (i.e. ,  A)  must  meet  condition  (46)  or 


tana 


4F2  -A2 
<  4FA 


(47) 


It  is  interesting  to  examine  the  case  of  a=0.  For  this  case,  the  image  of  the 
dish  is  a  circle  having  a  radius  of 


(y2  +Z2)1/2 
P  P 


4A 


F2  -FD 
4F2  -  A2 


As  is  expected  with  the  phased  array  at  the  focal  point,  that  is  F  =  D,  the  radius  is 
zero.  With  A  =  2F,  all  reflections  at  the  rim  are  perpendicular  to  the  dish  axis,  and 
the  image  is  infinite. 

The  path  length  from  the  dish  to  the  phased  array  is 

U  =  [<D-Qx)2  +  (yp-Qy)2  +  (zp“Qz)2)l/2  (48) 


The  positive  square  root  is  taken. 
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Ray  From  Phased  Array  to  Feed  Point 


Finally  the  phased  array  is  assumed  to  direct  the  reflected  ray  to  a  feed  point  on 
the  x  axis.  The  distance  to  the  feed  point  x=f  from  the  point  x=D,  y=yp,  z  =  Zp  is 


L,  =  [<D-f)2  +  y^  +  z^]1/2 


and  the  total  path  length  is  Eq.  (31),  (48),  (49) 

L  =  Li  +  L2  +  L3  (5 

Path  length  L  (or  for  that  matter  any  of  its  components)  can  be  converted  to  phase 
angle  by 


(p  =  — —  L  in  degrees 
A 


Where  X  is  the  wave  length  of  the  incident  radiation.  (X  and  L  must  be  in  the  same 
units.)  Equation  (51)  assumes  that  the  medium  in  the  antenna  system  is  homogeneous 
and  that  the  phase  changes  on  reflection  at  the  dish  and  at  the  phased  array  are  con¬ 
stants.  It  is  convenient  to  adjust  the  path  lengths  by  subtracting  one  L  from  the 
others;  since  only  relative  phase  is  required. 

This  treatment  of  the  antenna  system  does  not  consider  such  things  as  phase 
changes  at  reflection,  which  may  depend  on  the  angle  of  incidence;  polarization 
changes  at  various  places  in  the  antenna  system;  or  the  wave  nature  of  the  radiation 
that  leads  to  diffraction  effects  at  each  boundary.  These  phenomena  have  more  ef¬ 
fect  on  either  the  power  radiated  by  the  antenna  or  the  power  received  at  the  feed 
point.  However,  they  deserve  some  consideration  in  the  overall  understanding  of 
the  operation  of  the  antenna  described  here. 
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Appendix  A 


SUMMARY  OF  EQUATIONS  FOR  COMPUTATION 


The  numerical  computations  run  as  follows,  with  the  necessary  equations  re¬ 
peated  here  for  convenience: 


A.  Projections  of  Dish  and  Phased-Array  Rims  on  Wavefront 

The  coordinates  of  the  dish  rim  in  the  xyz  coordinate  system  are 


z  .  =  A  cos  9 

rim 

y  .  =  A  sin  0 

rim 


The  coordinates  of  the  projection  of  the  dish  rim  on  the  wavefront  plane  are 


y  .  cosa 
rim 


The  coordinates  of  the  phased-array  rim  are 


z' .  =  R„  cos  0 

rim  u 


y' .  =  R0  sin  0 

rim  0 


The  coordinates  of  the  projection  of  the  phased  array  are 


i 


It  is  sufficient  to  let  the  parameter  0  have  the  range 


-  x/2  <  0  <  x/2 

for  the  right  half  of  either  projection.  The  T  axis  is  taken  as  the  vertical  axis. 
Mirror  symmetry  about  the  Y  axis  gives  the  left  half  of  each  projection.  With  a  6 
increment  of  1*  or  2* ,  these  equations  give  sufficient  points  to  draw  smooth  curves 
with  the  automatic  plotters  used  with  digital  computers. 


B.  Image  of  the  Dish  Rim  on  the  Phased -Array  Plane 

Using  i  .  and  y  .  ,  the  ratios  of  the  components  of  the  unit  reflected  vector  R 
run  run  — 

at  the  dish  rim  are 


(2y* .  sin  a  -  4Fy  .  cos  a 

rim _ run _ 

(4F*-A*)  cos  a  -  4Fv  sin  a 


2(y  .  sina  -  2F  cos  a)  z  . 
rim _  rim 

(4F*-A2)  co s a  -  4F  y  .  sina 
nm 


(44) 


and  the  coordinates  of  the  intersection  of  the  reflected  rays  and  the  phased  -array 
plane  are 
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For  a  finite  rim  image  and  no  multiple  reflections,  the  following  conditions  must  be 
met 


R  >0 
x 


tana  < 


4F*  -A* 


Condition  (47)  is  for  a  finite  image,  and  Rx  >  0  is  for  no  multiple  reflection.  The  ray 
tracing  can  be  discontinued  at  this  point  if  either  of  these  conditions  is  not  met.  Again, 
an  increment  of  1*  or  2*  for  0  is  sufficient  to  draw  smooth  curves  for  the  dish  image 
on  the  phased -array  plane.  The  circle  which  is  the  phased  array  rim  can  also  be 
drawn  by  plotting  y^m  and  on  the  phased -array  plane. 


C.  Grid  of  Points  on  the  Phased  Array  Plane 


The  coordinates  of  the  grid  of  points  on  the  wavefront  plane  (which  are  the  start¬ 
ing  points  of  ray  tracing)  are  computed  from  the  following  equations.  These  points  all 
lie  on  or  within  the  projection  of  the  dish  rim.  The  number  of  Z  coordinates  is 


NZ  =  INT 


fc) 


This  is  the  number  of  points  for  one  half  of  the  Z  axis  including  the  origin  The  Z(J) 


Z(J)  =  £Z  (J-l);  J  =  l,2,  —  ,NZ 


The  number  of  Y(J)  coordinates  along  the  semiminor  axis  is 


NY(1)  =  INT 


r  A  cosot 


not  including  the  point  at  the  origin.  The  Y(T)  coordinates  are 
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Y(I)  =  AY '  (I-NY(l)-l)  - 


A7  sina 
4F 


(24) 


with 


I  =  1,2,  —  ,2NY(1)  +  1 


(25) 


For  Z(J)  (other  than  J  =  1) ,  the  range  of  the  index  I  is  found  from 


NY(J)  =  INT  cosa^ 


(22) 


and 


I  =  1  +NY(1)  -NYfcl),  —  ,1  +  NY(1)  +  NY(J)  (26) 

The  scheme  is  to  trace  a  set  of  rays  for  a  given  Z(J)  for  all  Y(I)  in  the  range  (26). 
Path  lengths  and  phases  for  this  set  of  Y(I)  can  then  be  listed  and  plotted  as  functions 
of  Y(I)  for  this  Z(J). 

D.  Tests  for  Rays  in  "Shadow  Area" 

Certain  of  the  Y(I)  will  fall  within  the  projection  of  the  phased  array.  If 

R,  s  ZP)  <  A  (16) 

then  all  Y(i)  with  indices  in  range  (26)  are  on  or  within  the  dish  rim  projection  and 
on  or  outside  of  the  phased -array  projection  ("shadow  area").  However,  if 

0  =£  Z(J)  <  Rj  (17) 


then  the  Y(I)  in  the  range 


cosa  -  D  sina  <  Y(I)  <  [Rj-Z2^)]1/2  -  D  sina  (18) 
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lie  within  the  shadow  area.  Rays  with  these  Y  (I)  can  still  be  traced  and  identified  so 
that  the  shadow  area  can  be  indicated  cm  the  various  plots  and  listings.  It  is  suffi¬ 
cient  to  note  the  smallest  and  greatest  I  meeting  conditions  (18)  for  each  Z(J). 


E.  Coordinates  of  the  Intersection  of  a  Ray  and  the  Dish 


Three  cases  are  considered 


(a)  a  =  0 


(b)  0  <  a  <  0.01  (radian) 

(c)  0.01  <  a  <  x/2  (radian) 


The  equations  for  Qy  are  (omitting  the  indices) 


(a)  Qy  =  Y,  a  =  0 


Qy  -  Y/cosa  -  |Z*  +  (Y/cosa)1! 


0  <  a  <  0.01 


cos  a 

(c)  Q  =  2F  . - 

y  sina 


r.  Z?sina+4FY  .  1*^1 

I1’!.1-  4F*co^a  SmaJ  ( 


0. 01  <  a  <  x/2 


For  all  a. 


QZ  =  Z 


Qx  = 


Q2  +Q2 

y  z 
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F.  Normal  to  Dish  at  Point  (Qx.Qy.Q^,  n 
The  components  of  the  unit  normal  are 

F 


nx  [FfF+Q^l1/1 


-Q. 


n  = 


y  2[F(F*Qx)l»/2 


(32) 


-Q. 


n  = 


z  2[F(F+Qx)]1/2 


G.  Unit  Vector  Along  Reflected  Ray,  R 
The  components  of  R  are 

R  =  (2H2  -1)  cos  a  t  2n  n  sin  a 
x  x  x  y 

R  =  (2nz  -1)  sina  +  2n  n  cos  a 

y  y  » y 

R  =  2n  fa  cosa  +  n  sir  a) 
z  z  x  y 

H.  Tests  for  Internal  Reflections 

(a)  A  necessary  condition  is  that 

R  >  0  (36) 

x 

If  this  test  is  not  met,  stop  ray  tracing  for  current  a,  make  note  of  this  fact,  and  go 
to  next  a. 
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(b)  With 


the  condition 


yd  *  R  ^/4F  ~V  *  **3 

X  J 


zd  >  f  (AV4F  -<y  ♦  Qz 


y2.  +  z  2  <  A2 
d  d 


must  be  met  for  the  reflected  ray  to  "escape"  from  the  dish  without  internal  reflec¬ 
tion.  If  it  is  not  met,  discontinue  ray  tracings  as  in  (a)  above. 


I.  Test  for  Ray  to  Intersect  Phased  Arra\ 


yP  -  r  (D-<y  *  % 


z  =  —  (D-Q  )  +  Q 
p  R  '  x  z 
x 


the  intersection  of  the  ray  and  phased-array  plane  lies  within  the  phased  array  if 


y2  +  z2  <  R§ 
P  F  0 


If  (42)  is  not  met,  ray  tracing  can  continue  but  the  ray  should  be  identified  as  being 
outside  the  phased  array.  After  all  rays  are  traced  for  a  given  a,  these  points  of 
intersection  can  be  plotted  on  the  phased-array  plans  along  with  the  image  of  the 
dish  and  the  rim  of  the  phased-array  plane  (see  Section  B). 
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J.  Path  Lengths  and  Phase 
The  path  lengths  are 


Ll  =  - 


<QX  cos  a  +  Q^sina) 


(31) 


h  =  l(D-Qj  +  (yp-Qy)2  +  (ZP^/11/2  (48) 

Ig  =  [(D-f)2  +  y*p  +  z^J1/2  (49) 

and  the  total  path  length  is 

L  =  L,  +  I,  +  L,  (50) 


This  path  can  be  adjusted  by  subtracting  some  constant  from  all  path  lengths.  The 
phase  is 


<P  = 


(51) 


provided  X  and  L  are  in  the  same  units.  If  L  is  in  feet  and  X  is  in  centimeters  then 

kfeet 

<p  =  10972. 8  — -  (degrees)  (51a) 

A 

cm 

Path  length  L  and  phase  <p  can  be  tabulated  and  plotted  as  functions  of  Y  for  each 
Z.  Note  which  rays  are  in  the  shadow  area  (see  Section  D). 
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Appendix  B 

COMPUTER  PROGRAM  OUTPUT 


The  ray  tracing  equations  developed  in  this  report  have  been  programmed  in 
FORTRAN  V  for  the  UNIVAC  1108.  Much  of  the  output  of  this  program  (named 
WAVE)  is  graphic  and  requires  the  Stromberg  DatagraphiX ,  Inc. ,  4060  plotter  and 
the  associated  INTEGRATED  GRAPHICS  SYSTEM  to  be  available  on  the  UNIVAC 
1108.  The  graphic  outputs  are  as  follows: 

First,  a  scale  drawing  of  the  antenna  system  showing  the  major  coordinate 
systems  is  made  with  the  plotter.  A  table  of  input  data,  such  as  the  focal  distance 
of  the  dish  and  the  radius  of  the  phased  array,  follows  this  drawing. 

Next,  a  plot  of  the  grid  of  points  on  the  wavefront  plane  is  made  along  with 
the  projection  of  the  rims  of  the  dish  and  the  phased  array.  Curves  of  phase 
change  as  functions  of  the  Y  coordinate  on  the  wavefront  plane  are  made  for  each 
Z  >  0.  Those  points,  which  lie  on  or  within  the  phased -array  rim,  are  marked 
with  the  symbol  *;  those  outside  are  plotted  with  the  symbol  • .  The  "shadow  area"; 
that  is,  those  rays  which  are  blocked  by  the  phased  array,  are  indicated  by  the 
shaded  area  on  the  plots.  On  another  plot,  the  points  of  intersection  of  the  rays 
with  the  phased -array  plane  are  drawn  cm  this  plane.  The  circle,  which  is  the 
phased-array  rim,  and  the  image  of  the  dish  rim  are  also  drawn.  In  this  plot, 
those  rays  in  the  shadow  area  are  not  shown. 

Printer  listings  of  the  path  lengths  and  phases  and  the  coordinates  of  the 
intersections  of  the  rays  with  the  phased-array  plane  can  be  made.  The  total 
number  of  rays  (those  in  the  shadow  and  those  on  or  within  the  phased  array)  can 
also  be  tabulated. 

Program  control  permits  selecting  the  several  outputs  individually  so  that 
only  those  which  are  required  can  be  produced.  Examples  of  graphic  outputs  are 
presented  on  the  following  pages.  The  charts,  which  are  exact  replicas  of  the 
computer  output,  have  been  reduced  to  85%  of  actual  size  to  conform  to  the  page 
size  of  this  report. 
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RADIUS  UK 

DISH  HIM, KT 

A  3 

•  .00 

KOCAL  DISTANCB.KT 

K  = 

4.00 

RADIUS  OK 

PHASED  ARRAY, KT 

RO  s 

1 .00 

LOCATION  OK  PHASBD  AKHAY , KT  D 
KBED  POINT, KT  KBBD 
*AVB  LBNOTM.CM  LAMBDA 


3. it 

0  oo 

5 . 31 


Computer -Generated  Scale  Drawing  of  Antenna  System.  The  dish  and  phased 
array  dimensions  and  locations  are  to  scale.  A  typical  wavefront  angle  a  is 
shown.  The  constants  in  the  table  hold  for  all  computer  plots  which  follow. 
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RADAR  RAY  TRAC  1 NU 

INPUT  DATA 


ANTENNA  CONSTANTS  IPT> 

RADIUS  OK  DISH  RIM  =  6.00 
KOCAL  DISTANCE  OK  DISH  =  4.00 
KEED  POINT  DISTANCE  =  .00 
DISTANCE  OK  PHASED  AREA  KROM  VERTEX  OK  DISH  :  3.50 
RADIUS  OK  PHASED  AREA  s  1.00 
LAMBDA  (WAVE  LENGTH  I  (CM  I  :  5.31 

NAVE  KRONT  PARAMETERS 

ALPHA  (INITIAL  ANOLE  OK  NAVE  KRONT  PLANE  s  .00 
INCREMENT  OK  ALPHA  (DEGREES!  :  3.50 
NUMBER  OK  ALPHAS  :  3 
INCREMENT  OK  Z  COORDINATE  ON  NAVE  KRONT  GRID(KT)  =  .40 
INCREMENT  OK  V  COORDINATE  ON  NAVE  KRONT  GRID(KT)  s  20 
PROJECTION  PARAMETER  INCREMENT  =  1.00 


Table  of  Constants  for  Antenna  Geometry  and 
Computer -Control  Parameters.  (A  plotter 
output) 
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-10.0  -1.5  -5.0  -2.5  -0.0  2.5  5.0  1.6  13-0 

AI.FMA  =  .00  COORD1NAT8  ON  WAVS  FRONT 

Is  1.00 


Phase  Angles  and  Path  Lengths  From  the  Wavefront  to  the  Feed  Point  as  Func¬ 
tions  of  the  Y  Coordinate  on  the  Wavefront.  This  plot  is  for  Z  =1.60  and 
a  =  0.  Note  that,  in  this  case,  no  rays  arc  in  the  "shadow  area".  The  symbols 
*  and  •  have  the  same  meanings  as  in  the  previous  plot. 


ALPHA  s  00 

A  -  5 ■ t‘ 


COORDINATE  ON  NAVB  PROMT 


Phase  Angles  and  Path  Lengths  From  the  Wavefront  to  the  Feed  Point  as 
Functions  of  the  Y  Coordinate  on  the  Wavefront.  This  plot  is  for  Z  =  5. 20 
and  =  0.  Note  that  all  rays  lie  outside  the  phased-array  rim  (plotted  with 
the  symbol  • ). 


AurriA  s 

*  1 


The  Phased- Array  Plane.  The  phased- array  rim  is  the  large 
circle.  The  smaller  one  is  the  dish  rim  image.  The  points  are 
the  intersections  of  the  rays  with  the  phased- array  plane.  The 
ray  intersections  in  the  "shadow  area"  are  not  plotted.  In  this 
case,  the  shadow  area  lies  in  a  circular  area  in  the  center  of  the 
phased- array  plane,  a  =  o. 


NWC  TP  4770 


rMIKTIOl  or  IIM 


ii  bins  on  un 


um  s  »  •• 


The  Wavefront  Plane,  the  Dish  and  Phased-Array  Rim  Projec¬ 
tions,  and  the  Grid  cf  Points  (+)  Used  to  Start  Ray  Tracing.  The 
smaller  ellipse  is  the  phased-array  rim  projection.  Rays  start¬ 
ing  on  or  within  this  ellipse  are  blocked  by  the  phased-array  and 
constitute  the  "shadow  area"  in  the  phase-angle  and  path-length 
plots.  a=  5.00. 
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ALPHA  =  5 . 00 
Z  :  »« 


COOftOI  MATE  ON  «AV8  FRONT 


Phase  Angles  and  Path  Lengths  of  Rays  From  the  Wavefront  to  the  Feed  Point 
as  Functions  of  the  Y  Coordinate  on  the  Wavefront.  This  plot  is  for  all  Y 
which  lie  along  the  line  Z  =  0. 40,  ot  =  5. 00  and  which  are  on  or  within  the  dish 
rim  projection.  The  shading  indicates  the  "shadow  area";  that  is,  the  region 
in  which  the  rays  are  blocked  by  the  phased  array.  The  symbol  *  indicates 
rays  which  are  reflected  so  that  they  intersect  the  phased-array  plane  on  or 
within  the  phased-array  rim.  The  symbol  •  indicates  rays  which  arrive  at  the 
phased-array  plane  but  lie  outside  the  rim. 
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ALTHA  s  5.55  COORDINATE  ON  RAVC  FRONT 

Z  =  Z.«0 


Phase  Angles  and  Path  Lengths  From  the  Wavefront  to  the  Feed  Point  as  Func¬ 
tions  of  the  Y  Coordinate  on  the  Wavefront.  This  plot  is  for  Z  =2.00  and 
a  =  5. 00.  Note  in  this  case  that  there  are  no  rays  in  the  shadow  area. 
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The  Phased- Array  Plane.  The  phased-array  rim  is  the  circle. 

The  image  of  the  dish  rim  is  the  cardioid  curve.  The  points  are 
intersections  of  the  rays  with  the  phased-array  plane.  The  "sha¬ 
dow  area"  is  seen  below  the  center  of  the  phased  array.  <*  =  5. 00. 
At  this  angle,  the  aberrations  of  the  antenna  system  become  appar¬ 
ent.  For  example,  in  the  region  near  YZ  origin,  there  is  no 
longer  a  unique  phase  for  a  given  point  on  the  phased-array  plane. 
In  other  words,  more  than  one  ray  from  the  wavefront  plane  can 
intersect  the  phased-array  plane  at  a  given  point  in  some  regions 
of  the  plane. 
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